arXiv:1503.04727vl [gr-qc] 16 Mar 2015 


Is motion under the conservative self-force in black hole spacetimes 
an integrable Hamiltonian system? 
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A point-like object moving in a background black hole spacetime experiences a gravitational self¬ 
force which can be expressed as a local function of the object’s instantaneous position and velocity, 
to linear order in the mass ratio. We consider the worldline dynamics defined by the conservative 
part of the local self-force, turning off the dissipative part, and we ask: Is that dynamical system a 
Hamiltonian system, and if so, is it integrable? 

In the Schwarzschild spacetime, we show that the system is Hamiltonian and integrable, to linear 
order in the mass ratio, for generic (but not necessarily all) stable bound orbits. There exist an energy 
and an angular momentum, being perturbed versions of their counterparts for geodesic motion, 
which are conserved under the forced motion. We also discuss difficulties associated with establishing 
analogous results in the Kerr spacetime. This result may be useful for future computational schemes, 
based on a local Hamiltonian description, for calculating the conservative self-force and its observable 
effects. It is also relevant to the assumption of the existence of a Hamiltonian for the conservative 
dynamics for generic orbits in the effective-one-body formalism, to linear order in the mass ratio, 
but to all orders in the post-Newtonian expansion. 


I. INTRODUCTION 

Recent years have seen great progress in computing 
the motion of and gravitational wave signals emitted by 
point-like objects orbiting massive black holes EHi , but 
some key foundational and computational challenges still 
lie ahead . A sufficiently accurate solution to the gravi¬ 
tational self-force problem will be crucial to the detection 
and analysis of signals from astrophysical extreme-mass- 
ratio inspirals by future space-based gravitational wave 
detectors [^. Study of the self-force also helps to inform 
other complementary approaches to the relativistic two- 
body problem, such as post-Newtonian theory and the 
effective-one-body formalism 

The foundations for computing the first-order self-force 
acting on a point mass in a curved background spacetime 
are provided by a result [3,S| which has been canonized as 
the ‘MiSaTaQuWa equation’. This is an equation for the 
self-force experienced by the point mass as a functional of 
its worldline history, arising from the tail part of its own 
linearized gravitational field. There are several different 
formulations of the MiSaTaQuWa result: 

1. A correction to the geodesic deviation equation: In 
this formulation, the deviation between the point 
mass’s actual accelerated worldline and a given 
fiducial background geodesic is characterized by a 
deviation vector field along the fiducial geodesic. 
The deviation vector obeys the geodesic deviation 
equation plus a self-force correction, and the self¬ 
force correction is constructed from the self-field 
sourced by the fiducial geodesic. This formulation 
is the version that is obtained directly from the 
most rigorous derivations 

2. A local equation of motion for an accelerated world¬ 
line: In this scheme, the worldline’s acceleration 
is given as a function of its instantaneous position 


and velocity, by using the geodesic which instan¬ 
taneously matches the position and velocity of the 
accelerated worldline (the osculating geodesic) as 
the source to compute the instantaneous self-force 

i- 

3. A non-local equation of motion for an accelerated 
worldline: In this family of schemes, one uses the 
history of the actual accelerated worldline as the 
source for the self-force^. Viewing the self-force 
as a given functional of the worldline, this corre¬ 
sponds to a non-local equation of motion for the 
the worldline—the instantaneous acceleration de¬ 
pends on the worldline’s entire history. 

Schemes 2 and 3 have a clear advantage over scheme 
I in being able to describe inspiraling motion over a ra¬ 
diation reaction timescale. Both 2 and 3 reduce to 1 for 
small deviations from the fiducial geodesic, over short 
timescales. Scheme 2 is said to arise from 3 via a reduc¬ 
tion of order procedure: because the self-forced worldline 
and the osculating geodesic agree to zeroth order in the 
mass ratio, the former can be replaced with the latter in 
calculating the self-force to first order. The relative accu¬ 
racies and computational merits of 2 and variants of 3 are 
subtle issues, entangled with formulating and implement¬ 
ing a consistent second-order analog of the MiSaTaQuWa 
result (lll - [T3| . which will not be addressed here. 

This paper addresses a formal question concerning the 
worldline dynamics defined by scheme 2, using only the 
conservative part of the self-force, turning off the dissipa¬ 
tive/radiative part. We ask: Is the resultant dynamical 


^ Since the stress energy tensor for an accelerated worldline is not 
conserved, implementing this sheme requires an ad hoc relaxation 
of the Lorentz gauge condition. 
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system, for the worldline degrees of freedom only, an in- 
tegrable Hamiltonian system, to linear order in the mass 
ratio? 

In the case where the background spacetime is 
Schwarzschild, we show that this system is indeed a 
Hamiltonian dynamical system, in the classic sense [l5l| . 
We also show that it is integrable; in other words, there 
exist three functions (the rest mass, an energy, and an 
angular momentum) on the effectively six-dimensional 
phase space whose values are instantaneously conserved 
along the system’s trajectories. These results hold to 
linear order in the mass ratio, for generic stable bound 
orbits (with some caveats for orbits in the zoom-whirl 
regime). 

We also carry out a partial analysis of the case of 
the Kerr spacetime, following the same line of reasoning 
which leads to a successful construction in Schwarzschild. 
In Kerr, we encounter complications which are associ¬ 
ated with orbital resonances [l^, when the ratio of the 
frequencies of the radial and polar motions is a rational 
number. These complications prevent us from drawing 
definite conclusions about the Kerr case. 


There is, however, a heuristic argument which suggests 
that the motion in Kerr should be integrable and Hamil¬ 
tonian [l3, [13 ■ Namely, the ambiguities in the defini¬ 
tion of angular momentum related to the BMS group are 
associated with the dissipative sector of the linear per¬ 
turbations, so the conservative sector should admit three 
independent conserved components of angular momen¬ 
tum [H, a sufficient number to make the system inte¬ 
grable. This general picture is consistent with what has 
been found in the post-Newtonian approximation, where 
the system is integrable at successive orders [18 1. 

The conservative self-force dynamics in Kerr has re- 
centl y b een cast in a Hamiltonian-like formulation in 
Ref. [^. However, that formulation is not Hamiltonian 
in the classic sense used here. Specifically, in that formal¬ 
ism, the equations of motion are obtained by differentiat¬ 
ing the Hamiltonian with respect to one set of variables 
while holding a second set of variables fixed, and then 
taking the two sets of variables to coincide. 

The existence of a Hamiltonian system for the con¬ 
servative dynamics is a foundational assumption in the 
effective-one-body formalism [2l|, as well as in other re¬ 
lated treatments of the relativistic two-body problem, 
e.g. [22l - [^ . While the results of Refs. [23i - l^ make it 
clear that such a Hamiltonian exists for circular orbits in 
the extreme-mass-ratio limit (to linear order in the mass 
ratio, but to all orders in the post-Newtonian expansion), 
our result confirms the validity of this assumption for 
generic (non-circular) orbits, in the non-spinning case. 


This paper is organized as follows. Our construction 
relies heavily on the fact that geodesic motion in the 
Schwarzschild (and Kerr) spacetime is completely inte¬ 
grable, and thus admits a representation in terms of gen¬ 
eralized action-angle variables [l3, . We review rele¬ 

vant properties of these variables in Sec. El and we review 


relevant properties of the local conservative self-force in 
SeaUm We develop sufficient conditions for the forced 
system to be Hamiltonian and integrable in Sec. IIVI and 
we conclude in Sec. El 


II. GEODESIC MOTION IN KERR AND 
ACTION-ANGLE VARIABLES 

Geodesic motion in a spacetime with metric g^i, is gen¬ 
erated by the Hamiltonian 

( 2 . 1 ) 

where z^{X) are the coordinates of a worldline, and Pfj.{X) 
are the components of its momentum, together with the 
canonical symplectic form Hn = dz^ A dp^ on the world¬ 
line phase space {z^,Pfi) [la[i3- Hamilton’s equations 
read 

^ ^ dpij, ^ dHp 

dX dp^ ’ dX dz^^ ’ 

and imply the geodesic equation in affine parameteriza¬ 
tion, p'^VuPii = 0. Identifying the affine parameter as 
A = t/to, where r is the proper time along the world¬ 
line, gives the usual expression p^ = mgfj^^dz'^/dr for 
the momentum of a particle of mass m, and also yields 
Hq = —to^/ 2 on shell. 

In the Kerr spacetime, with Boyer-Lindquist coordi¬ 
nates z^ = geodesic motion is completely in¬ 

tegrable thanks to the existence of four first integrals of 
motion. These are the Hamiltonian i7o, the energy E = 
— {dt)^P)j., the axial angular momentum = {d^)^Pfj,, 
and the Carter constant Q = K^'^Pf^p,,, where (dt)^ and 
are the timelike and axial Killing vectors, and 
is the non-trivial Killing tensor. 

These four first integrals are independent and in invo¬ 
lution, which implies the existence of generalized action- 
angle variables^ for bound geodesics in the Kerr geome¬ 
try [ 13 , [2^. These are canonical coordinates (g“, Jq) = 
(g*, q^, , q'^, Jt, Jr, Je, J<j>) on the worldline phase space, 

for which the geodesic Hamiltonian Hq depends only on 
the action variables Ja, 

Hoiz,p) = HoiJ), (2.3) 

and not on the angle variables g“. They are obtained 
from the {z^,p^) coordinates via a canonical transfor¬ 
mation 

= Jr,{zf^,pf,), (2.4) 


^ The action angle variables discussed here are associated with 
the affine parameter X = rim along the worldline; there are also 
other action-angle coordinates on the phase space associated with 
Mino time [5^ and Boyer-Lindquist coordinate time. 
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so that r^o = = dgaAdJc,. Hamilton’s equations 

then take the particularly simple form 


IX 

dJo^ 

~dX 


OHq _ 
dJc " 
dH^ 

dq°‘ 


= 0 . 


(2.5a) 

(2.5b) 


The action variables Ja are all independent constants of 
motion, and the angle variables g“ all increase linearly, 
at the constant rates known as the fundamental fre¬ 
quencies. The angle variables , q^, and q’^ are each pe¬ 
riodic with period 27r, while g* has an infinite range. The 
action variables Ja are functions of the geodesic first in¬ 
tegrals Pa = {Ho,E,L 2 ,Q)', in particular, Jt = —E and 

Jcl> — L z • 

In the Schwarzschild limit of the Kerr geometry, both 
geodesic motion and self-forced motion are confined to a 
plane, which can be taken without loss of generality to 
be the equatorial plane 9 = 7r/2. We can then ignore 
the 0-motion, working in the reduced phase space with 
coordinates (z^,p^) = We have the 

three first integrals Pa = {Ho,E,Lz), and action-angle 
variables (g“, Ja) = (g*, g’’, q'^, Jt, Jr, J<t,), defined just as 
in the general Kerr case. 


III. CONSERVATIVE-SELF-FORCE 
PERTURBATION TO GEODESIC MOTION 


and the four variables Ja [13: 

r{q,J) = (3.3a) 

Pa{q,J) = Fa{q\q\j). (3.3b) 

In the Schwarzschild case, these functions depend only 
on q'" and the three Ja- 

The forcing functions /“ and Fa are periodic functions 
of q^ and of g®, each with period 27r. They can thus be 
expanded as Fourier series in q’’ and g®, according to 

OO 

Fa{q\q\j)= Y. Fa{K,ke,J)E^^‘'^+^^o^\ 

kr ,kQ — — OC 

(3.4) 

and similarly for /“. As shown by Mino [^, the (0,0) 
Fourier mode of each Fa vanishes, 

1 p27T p27r 

Fai0,0,J) = -^ dq^ dg®F„(g^g^J) = 0, 
(Stt)^ Jo Jo 

(3.5) 

due to reflection properties of Kerr geodesics and to the 
time-reversal symmetry of the conservative self-force [T^ . 

|23. 

In the Schwarzschild case, we lose the dependence on 
the 0-motion. Equations (13.31) and (13.41) are replaced by 

OO 

Fa{q\J)= Y FXkr,J)F’‘^‘>\ (3.6) 


Instead of geodesic motion, we now consider a point 
mass m with worldline z^(A) experiencing a local linear 
perturbing force, 

(j 

-_=pA‘, p’'y^p^ = ePij,{z,p), (3.1) 

where e is a small parameter. We are interested in the 
case where the forcing function P^{z,p) is given by the 
conservative part of the osculating-geodesic-sourced first- 
order gravitational self-force in the Kerr spacetime, with 
the parameter e being the small mass ratio m/M. 

The worldline dynamics can be re-expressed in 
terms of the geodesic action-angle variables ((/“, Ja) as 


(jn^ 

^ = oj^{J)+er{q,J), (3.2a) 

^ = eFa{q,J), (3.2b) 


following Ref. [17[. Here we use the same phase space 
coordinate transformation (1^ as for geodesic motion 
to obtain Eqs. dSSI) from Eqs. (I3T]). The forcing func¬ 
tions /“(g, J) and Fa{q, J) are determined from the self¬ 
force components P^{z,p) via /“ = (9g“/9p^)^ and 
Fa = idJa/dPfi)z Ffj,. These functions have the impor¬ 
tant property that they are independent of the angle vari¬ 
ables g* and q'^, because of the symmetries of the Kerr 
spacetime. Thus, they can written as functions of q^, q^ 


and similarly with Fa replaced by /“. Equation (j3.5l) 
becomes 


FX0,J) = ^ [ dq^ Fa{q'^,J) = 0. (3.7) 

27r Jo 

Recalling that the forcing functions Fa give the rates 
of change of the action variables Ja, Eqs. (1^ and (Irfl) 
express the fact that the conservative first-order self-force 
causes no net change in the geodesic first integrals, when 
the force is evaluated along a geodesic and suitably aver¬ 
aged [ 13 . In the Schwarzschild case (1X71) . the averaging 
is an orbital average or time average over one period of 
radial motion. In the Kerr case (13.51) . the average is over 
the (g’",g®) torus in phase space, which is equivalent to 
a time average only over an infinite time and only for 
non-resonant orbits [H, [T3, [13 ■ 


IV. IS THE PERTURBED SYSTEM 
HAMILTONIAN AND INTEGRABLE? 


The perturbed system (1^ will be Hamiltonian and 
integrable, to linear order in e, if there exist new phase 
space coordinates (g“, Ja) and a new Hamiltonian func¬ 
tion H{J) for which Eqs. (13.21) are equivalent to 




dH{J) 

dJa 


00 ’). 


§ = (4.1) 
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Without loss of generality, we can express the new coordi¬ 
nates as linear perturbations of the geodesic action-angle 
coordinates (g“, Jc): 

r{q,J) = g“ + ex“(g,^), (4.2a) 

Ja{q,J) = Ja+ eCa{q,J), (4.2b) 

for some functions x“ and Cq, to be determined. Note that 
g3) is not assumed to be a canonical transformation. 
Similarly we can express the new Hamiltonian as 

H{J) = HoiJ) + eHiiJ), (4.3) 


where Hq is the geodesic Hamiltonian function, for some 
function Hi{J) to be determined. 

Combining Eqs. (14.11) . (14.21) and (14.31) now yields that 
Eqs. (14.11) will be equivalent to Eqs. (13.21) to linear order 
in e if 


r{q,J) 

Fa{q,J) 


dCa 

dq!^' 


duj°‘ dHi{J) 


(4.4a) 


(4.4b) 


where a;“ = a;“(J) are the geodesic fundamental fre¬ 
quencies (1231) . Thus, the conservative-self-force dynam¬ 
ics (133) will be Hamiltonian and integrable if there 
exist functions x“(g, J), Ca{q,J), and Hi{J) satisfying 
Eqs. (I4.4al) and (I4.4bl) . 

In light of Eqs. (13.31) . it is natural to consider solutions 
for x“ and (a which, like /“ and Fa, are independent of 
g* and q^, and which are periodic functions of g’' and g® 
[or of just g’’ in Schwarzschild]. We can then decompose 
all of these functions into discrete Eourier series for the 
g’' and g® dependence, just as for Fa in Eq. (EH) [or 
Eq. (13.61) ]. This defines Eourier mode amplitudes /“, Fa, 
X°‘, and Ca which are functions of the two integers kr and 
ke [or just kr] and all of the action variables Jq,. 

We then have the following Fourier transforms of 
Eqs. (I4.4al) and (I4.4bl) : 


f)l ,OL f)TJ 

r = -*(a;-A:)r+ ^C/3 + 4..o4.,o^, (4.5a) 
Fa = -i(uj . k)(a, (4.5b) 


where 
(w • k) 


uj'^kr + Lo^ke 
oj'^kr 


Kerr 

Schwarzschild ’ 


(4.6) 


and with Skgfi ^ 1 in Schwarzschild. If we restrict at¬ 
tention to Fourier modes for which {uj ■ k) ^ 0, then 
Eqs. (I4.5al) and (|4.5bl) admit the simple solutions 



iFg 

{uj ■ k)’ 


X 




(4.7) 


In the general Kerr case, the quantity uj ■ k = kr + 
uj^kg can vanish at locations in phase space where juj^ 
is a rational number, corresponding to an orbital reso¬ 
nance in the r and 9 motions [l^ . The solutions ( 1371 ) 


are clearly not valid in such cases, and so our analysis 
does not allow us to draw any definite conclusions about 
the Kerr case. 

For stable bound orbits in Schwarzschild, the quan¬ 
tity (jj ■ k = u}''kr vanishes only when kr = 0, since 
w’' = 0 occurs only in the limit of unbound or unsta¬ 
ble orbits. Equations ( 1371 ) thus provide valid solutions 
for all Fourier modes of and (a, except for the kr = 0 
modes. Given the fact ( 1371 ) that Fa{0,J) = 0, we see 
that a separate solution to Eqs. (14.Sal) and (14.5bl) for the 
case kr = 0 is given by 

Co(0,J)=(^) F{0,J), H,=0, (4.8) 

[with x“(0, J) unconstrained], provided that the matrix 
doj^ fdJa is invertible. It follows from the results of 
Ref. that duj^/dJa is invertible for all stable bound 
geodesics, except for those along the singular curve 
associated with isofrequency pairing of Schwarzschild 
geodesics in the zoom-whirl regime (and those along the 
separatrix defining the boundary of stable orbits); see 
Figure I of Ref. [^. Thus, we have constructed a so¬ 
lution of Eqs. (I4.4al) and (14.4bl) . and so the perturbed 
motion is Hamiltonian and integrable. 


V. CONCLUSION 

We have shown that the local first-order conservative 
self-force dynamics in the Schwarzschild spacetime is an 
integrable Hamiltonian system, to linear order in the 
mass ratio, for generic stable bound orbits outside the 
zoom-whirl regime (more specifically, for all orbits to the 
right of the “singular curve” in Figure 1 of Ref. [1^). The 
Hamiltonian system is defined by Eqs. EH, with the co¬ 
ordinates (g“, Jq) defined in terms of the geodesic action- 
angle coordinates by Eqs. ( 133 . and with the Fourier 
modes of the functions and Ca given by Eqs. (133 
and (14.81) . The quantities —Jt and J^j, are well-defined 
functions on the worldline phase space, which are per¬ 
turbed versions of the geodesic energy E and angular 
momentum L^, which are instantaneously conserved to 
linear order under the conservative-self-forced motion. 

Finally, we remark that the obstacles encountered by 
our construction in the general Kerr case do not show 
that the system is not Hamiltonian or integrable in Kerr. 
As mentioned in the introduction, there is a heuristic 
argument indicating that the dynamics should be inte¬ 
grable and conservative in Kerr. 
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